Soliton dynamics in deformable nonlinear lattices 
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We describe wave propagation and soliton localization in photonic lattices which are induced 
in a nonlinear medium by an optical interference pattern, taking into account the inherent lattice 
deformations at the soliton location. We obtain exact analytical solutions and identify the key factors 
defining soliton mobility, including the effects of gap merging and lattice imbalance, underlying the 
differences with discrete and gap solitons in conventional photonic structures. 
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The effect of Bragg scattering from a periodic poten- 
tial is a fundamental phenomenon which is responsible 
for a strong modification of wave dispersion and the ap- 
pearance of spectral gaps. The structure of band-gap 
spectrum in crystals defines the electron transport prop- 
erties, and similar concepts were successfully developed 
in the field of optics . The ultimate flexibility in man- 
aging wave transport and localization may be achieved 
in dynamically induced lattices. Such reconfigurable lat- 
tices can be realized in any nonlinear media, where a 
modulated wave can modify the medium characteristics 
(e.g. an optical refractive index E3- E3| ) and induce an 
effective periodic potential. 

Nonlinearity also supports wave localization inside the 
band-gaps in the form of discrete and gap solitons, 
when dispersion or diffraction is suppressed through self- 
focusing H, IE S S @- The dynamics of such soli- 
tons in fixed lattices is strongly influenced by a self- 
induced Peierls-Nabarro (PN) potential which can result 
in soliton trapping, and this effect has applications for 
intensity-dependent beam steering |lfjj . Recent theoret- 
ical studies have indicated that the PN potential can 
vanish in nonlinear lattices which become deformed at 
the soliton location 0, 0], however soliton trapping 
was observed in recent experiments due to soliton-lattice 
interaction Q. In this Letter, we describe the key 
mechanisms which determine the soliton mobility in de- 
formable nonlinear lattices when the PN potential is ab- 
sent. In particular, we identify the fundamental effect of 
band-gap merging, and suggest how it can be controlled 
through the lattice imbalance. 

We will analyze the case when the interaction between 
the localized beam and the nonlinear lattice is phase- 
insensitive. This situation is realized for optical waves 
which are mutually incoherent 0, LH Ell ■ We consider 
the (1+1)D geometry, where waves can diffract in one 
spatial dimension, and they are confined by a guiding 
potential in the other transverse dimension. Then, the 
wave dynamics can be approximately described by the 
coupled nonlinear Schrodinger (NLS) equations for the 
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where x is the transverse spatial coordinate, z is the prop- 
agation direction, / = |i?i| 2 + l-E^I 2 is the total inten- 
sity, and a = ±1 stands for the focusing or defocusing 
nonlinearity, respectively. We have neglected the higher- 
order nonlinear effects in Eq. (|TJ to identify the most 
general ph ysic al effects. The model Eqs. Q are fully 
integrable [14j. and their solutions can be obtained ana- 
lytically. 

We first summarize the properties of nonlinear lattices 
which are found as stationary periodic solutions of 
Eq. (JTJ in the most general form, E\ = ri(x)e lVl ^ x ' +ll3lZ . 
Here the amplitude [r(x)] and phase [<£>(x)] profiles are 
fixed, and the propagation constant /3i is proportional 
to the wavevector component along the z direction. We 
can present the lattice intensity profile in the same form 
for both the cases of self- focusing (a = +1) and self- 
defocusing (<r = —1) nonlinearities, 



ar 2 (x) — A 2 cn 2 (xn, m) + Vq, 
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where cn is the Jacobi elliptic function with modulus m 
(0 < m < 1), F = Pi/3 - K 2 (2m ~ l)/3, A = Ky/m, 
K = 4K (\/m) j (2d) , K is the complete elliptic integral 
of the first kind, and d is the lattice period. The lattice 
phase is given in an integral form, tpi (x) — J Cr^ 2 (x) dx, 
where C 2 = -V (V + A 2 )(V + A 2 - k 2 ). The lattice is 
defined by the following parameters: period d, modulus 
m, and the propagation constant (3\, which should be 
chosen to satisfy the conditions A 2 > and C 2 > 0. 
The nonlinear waves exhibit strong modulational insta- 
bility when a = +1 and Vq ~ k 2 (1 — to). In the fol- 
lowing we consider the lattices with parameters a = +1, 
Vo — and a = — 1, Vq ~ — tor 2 , when only weak os- 
cillatory instabilities may appear, which are suppressed 
in media with saturable nonlinearity [ill . Such so- 
lutions can be excited experimentally by two interfer- 
ing waves representing the dominant Fourier harmonics, 
E\(x,z = 0) ~ F+ exp(iTTx/d) + F- exp(— iirx/d). The 
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FIG. 1: Top rows: Characteristic profiles of nonlinear lattices 
with d = 5, m = 0.3, and the same profiles err? (a;) (top), but 
different nonlinearities and phase structures: (a) flat-phase in 
a self-focusing medium (C = 0, Vo = 0, a = +1), (b) flat- 
phase in a self-defocusing medium (C = 0, Vo = — ^4 2 , o = 
— 1), and (c) nontrivial phase modulation in a self-defocusing 
medium (C ~ 0.12, V = -^ 2 -0.1, a = -1). Bottom: Linear 
Bloch-wave dispersion, the same for all the lattices (a-c). 



balanced excitation {F_ = F+) produces a lattice with 
C = and a flat-phase profile [see Figs. [IJa,b)], and im- 
balance leads to nontrivial phase modulation with C =/= 
[Fig. He)]. 

The periodic modulation induced by a nonlinear wave 
strongly affects dynamics of the probe beam which can 
exhibit Bragg scattering from the lattice. Since the 
beam-lattice coupling is phase-insensitive, the propaga- 
tion of small-amplitude probe beam (with an envelope 
£2) is governed by linear equation with a stationary pe- 
riodic potential defined by the unperturbed lattice inten- 
sity profile [Eq. J3J)] , and therefore does not depend on the 
lattice imbalance. Linear wave dynamics in a periodic po- 
tential can be analyzed by decomposing the wave packet 
into a superposition of extended eigenmodes called Bloch 
waves which are found as solutions of the wave equation 
in the form E%{x,z) = B(x) exp(ij%z + K^x/d), where 
K is the normalized Bloch wave number and B(x] = 
B(x+d) is a periodic wave profile. We immediately notice 
that nonlinear periodic waves [E\(x, z)] satisfy the Bloch 
condition. Moreover, for fixed parameters d and m, the 
lattice intensity profiles are exactly the same [up to a con- 
stant shift defined by the value of Vo(/?i), see Fig. a-c)], 
and these solutions form a full set of Bloch waves which 
wave- numbers are found as K — <p(d) — <p(0). Conversely, 
for all lattices with particular d and m but different phase 
structure and type of nonlinearity, the dispersion of Bloch 
waves Pb{K) is equivalent, see Fig. ^bottom). There 
exists a Bragg-reflection gap for a range of propagation 
constants n 2 {m - 1) + 2V < (3 b < K 2 (2m - 1) + 2V , 



and a semi-infinite gap due to total internal reflection for 
ft > K 2 m + 2V . 

In order to uncover the fundamental features of non- 
linear wave transport in deformable lattices, we analyze 
the properties of stationary and moving solitons, con- 
sisting of nonlinearly coupled spatially localized and ex- 
tended lattice components. We consider the most general 
case of phase-modulated lattices with internal imbalance 
(C^0), and identify new effects in comparison with soli- 
tons in flat-phase lattices [ill fl2T | . We use the approach 
based on Darboux transformation to "add" a local- 
ized soliton to a periodic lattice. In this method, first 
one calculates auxiliary functions as solutions of equa- 
tions associated with the Lax pair representation of the 
integrable Eqs. QJ, 
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Here £ is the complex parameter which implicitly de- 
fines the soliton amplitude, width, and speed of its mo- 
tion across the lattice. The functions ip n can be found 
analytically for lattices with flat phase (C = 0) pro- 
files I n the general case of lattices with arbi- 
trary phase structure (C 7^ 0), we can still derive a 
number of key analytical relations. We note that the 
equations for ipi and 02 in Q are uncoupled from 
03 and their general solution can be represented as a 
sum of two eigenmodes. With no loss of generality, we 
choose the eigenmodes which profiles remain stationary 
along the z direction, similar to the underlying lattice, 

4>il(x, z) — ipil( x ) ex P(i7i 2*0) where j — 1,2 is the 
eigenmode index. The propagation constants are found 
as 7 « = 7 « + ft = ft/2 ± [(-ft + 2a\E 1 \ 2 + £ 2 ) 2 + 

4a(dEi/dx + £Ei) {dE{/dx - ^*)] 1/2 , and 73 = -£ 2 / 2 - 
There is a specific relation between the amplitudes of 
two components, ip[ 3) /ip { 2 3) = (dEi/dx + ^i)/(~7^ ) + 
£|£a| 2 + £ 2 /2). Using these constraints, we can find the 
profiles of eigenmodes by integrating Eqs. J3J numeri- 
cally. Then, the solution of original Eqs. (|TJ describing a 
composite soliton on a lattice is found as 



Ei(x,z) = E 1 - 2Re(f)0i0^L>- 1 , 
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where E± is the nonlinear lattice profile, and D = 
ffl^il 2 + 1 02 1 2 + 103 1 2 - We can reveal the physical ori- 
gin of solutions associated with different eigenmodes tp^' 
by explicitly taking into account soliton evolution along 
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the propagation direction z, and write its profile in the 
following equivalent form, 
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where x( x ) = m V <7 IV , i" 
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(5) 
and 



the soliton propagation constant is p$ — ^\ J> — 73. It 
follows from Eq. J^J that the average velocity of soliton 
motion across the lattice is u — lim( a;2 _ Xl )_ +00 Im(73 — 
li){x2 — Xx)/[x{x2) — x{ x i)]- The initial soliton position 
xq at z = (or a fixed location of the stationary soliton) 
can be found from the condition x( x o) — ; and it de- 
pends on the value of tps, which is an arbitrary constant 
according to the last two equations in ©• Since the soli- 
ton is localized, 02 should belong to a band-gap, where 
linear waves cannot propagate and become trapped at the 
self-induced nonlinear defect. We find that for stationary 
solitons (v = 0), the propagation constants (3% associ- 
ated with two eigenmodes (j — 1,2) of Eqs. J2J appear 
inside the total internal reflection or the Bragg-reflection 
gap, respectively, and the specific value of 02 depends on 
the soliton parameter £. Accordingly, solitons can exist in 
both of these gaps in case of self-focusing (er = +1) non- 
linearity. On the other hand, only Bragg-gap solitons can 
form in media with self-defocusing (er = —1) nonlinear- 
ity, as the denominator D in Eq. Q becomes singular for 
modes in the other gap. These conclusions are in agree- 
ment with earlier numerical results for flat-phase non- 
linear lattices As a matter of fact, these existence 
properties may seem to be the same as for solitons in 
fixed lattices (see Ref. ^j) , and references therein) , how- 
ever there is a fundamental difference. For every value 
of the propagation constant, the position of stationary 
solitons with respect to the nonlinear lattice can be ar- 
bitrary depending on the choice of ipg [cf. two bottom 
rows in Figs. Eta, c) and Eta)], whereas only two specific 
locations on a period (one stable and one unstable) are 
possible in fixed lattices 0] ■ This indicates the absence 
of the self-induced Peierls-Nabarro potential-which can 
inhibit soliton motion through fixed lattices [1(1 ]. 

Even in the absence of the Peierls-Nabarro potential, 
certain conditions have to be satisfied to sustain soliton 
motion through a nonlinear lattice. We note that the 
lattice is deformed at the soliton location, and if the cor- 
responding defect is not created at the input the soli- 
ton mobility may be restricted. Even for optimal in- 
put conditions, we find that there appear fundamental 
limitations on the structure and mobility of solitons re- 
lated to the properties of the band-gap spectrum. We 
scan the full complex plane of soliton parameters £ and 
calculate how the existence regions of soliton propaga- 
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FIG. 2: Top: Existence regions of the existence of station- 
ary and moving composite band-gap solitons (shaded) in a 
self-focusing medium (a — +1); the normalized propagation 
constant is @2 — 2Vo- Bottom rows: (a-c) Intensity profiles 
of the bright soliton (solid) and lattice (dashed and shaded) 
components corresponding to marked points in the top plot. 
Two rows illustrate different transverse positions. Lattice pa- 
rameters correspond to Fig. 0a). 



tion constants {(3 2 ) change with the variation of their 
velocity. We first consider the properties of solitons in 
flat-phase lattices, and then use these results to iden- 
tify phenomena due to lattice phase modulation. In the 
self- focusing case (er = +1), the existence regions expand 
at large velocities [see Fig. Et top)]. We obtain that the 
moving solitons can exist everywhere inside the dynam- 
ical bandgaps of the lattice spectrum, which edges are 
found from the linear Bloch-wave dispersion relations as 
(—ddf3bl dKb, (3b). The dynamical gap broadening is a 
fundamental phenomenon, however earlier studies of this 
effect were based on a simple coupled-mode theory which 
accounts for a single isolated gap (see Ref. 0], and ref- 
erences therein). Our results describe a physical system 
with a non-trivial multi-gap spectrum, and we uncover 
the key effect of gap merging at large soliton velocities. 
In this regime, the propagation constant is detuned from 
the resonance with the lattice, and the large-amplitude 
soliton begins to locally erase the lattice, see Fig.|2tb). 

The features of moving bright solitons in media with 
self-defocusing nonlinearity are fundamentally different 
from the self-focusing case. At relatively small ve- 
locities, we observe the expansion of existence region 
which fully occupies the dynamic Bragg-reflection gap 
[see Fig. [3], however at larger velocities the existence re- 
gion shrinks and completely disappears above a critical 
velocity. This happens because the soliton localization 
in self-defocusing media is only possible in the Bragg- 
reflection gap, which disappears through gap merging, 
as the effectiveness of scattering is reduced at larger ve- 
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FIG. 3: Soliton existence region and profiles in a self- 
defocusing medium (a — — 1). Notations are the same as 
in Fig. |21 and lattice parameters correspond to Fig. 0^b) . 
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FIG. 4: Top: Soliton existence region for a lattice with 
nontrivial phase modulation [Fig. ^c)] in a self-defocusing 
medium shown with shading. Dashed line: boundary of the 
existence region of a trivial-phase lattice [Fig. Q^b)]. Bottom 
rows: Intensity profiles of the lattice and soliton components 
corresponding to the marked points in the top plot. 



locities. Accordingly, only small-amplitude solitons are 
supported close to the existence boundary [see example 
in Fig. He)]. 

We now analyze the effect of lattice imbalance and as- 
sociated phase modulation on the soliton dynamics. We 
have established that imbalanced lattices still have sym- 
metric intensity profiles [see Eq. 10 and Fig.^c)]. Since 
the coupling of spatially localized and periodic compo- 
nents is phase-insensitive, the left and right propaga- 
tion directions are exactly equivalent for linear (small- 
amplitude) wave packets in the second component. How- 
ever, at larger wave amplitudes, i.e. in the regime of soli- 
ton formation, the lattice becomes locally deformed and 
this process does depend on its phase, resulting in strong 
differences between the left- and right- moving solitons 



[cf. Figs. E{b) and (c)]. This is clearly indicated by a 
dramatic change in the soliton existence region in a self- 
defocusing medium [Fig.^top)], which becomes strongly 
asymmetric. For comparison, in the same plot we show 
with the dashed line the boundary of the existence re- 
gion for a flat-phase lattice. We see that the existence 
regions exactly coincide and cover the whole dynamic 
Bragg-reflection gap up to a critical velocity. At large 
velocities, the gap-merging effect comes into play and 
existence region critically depends on the lattice phase: 
solitons can travel at much large velocities in the direc- 
tion defined by the lattice phase gradient. We stress that 
this is a nontrivial result, which has no analogs for fixed 
periodic structures. 

In conclusion, we have described the dynamics of soli- 
tons in deformable lattices which are induced by periodic 
waves in a Kerr-type nonlinear medium. We have used 
exact analytical solutions to demonstrate that, even in 
the absence of the Peierls-Nabarro potential, the soliton 
mobility can be restricted due to other physical mecha- 
nisms, including the effect of band-gap merging. We have 
also suggested a new approach for controlling soliton mo- 
tion by introducing a lattice excitation imbalance, while 
at the same time preserving the linear wave spectrum. 

A.S. thanks Yuri Kivshar, Anton Desyatnikov, and 
Evgueny Doktorov for useful discussions and comments. 
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